arXiv: 1502.05285v2 [cond-mat.mtrl-sci] 22Jun2015 


Structure and dynamics of the fullerene polymer L^Ceo 
studied with neutron scattering 

S. RolsA^D. Pontiroli, 2 , [j] C. Cavallari , 1 M. Gaboardi , 2 M. Aramini , 2 D. Richard , 1 
M. R. Johnson , 1 J. M. Zanotti , 3 E. Suard , 1 M. Maccarini, 1 ’^ and M. Ricco 2 

1 Institut Laue Langevin, BP 156, 6 rue Jules Horowitz, 3804-2 Grenoble Cedex, France 
2 Dipartimento di Fisica, Universita di Parma, 

Via G. Usberti 1/a, 43100 Parma, Italy 
3 Laboratoire Leon Brillouin, CEA Saclay, 

91191 Gif Sur Yvette Grenoble Cedex, France 
(Dated: June 23, 2015) 

Abstract 

The two-dimensional polymer structure and lattice dynamics of the superionic conductor LRCgo 
are investigated by neutron diffraction and spectroscopy. The peculiar bonding architecture of 
this compound is confirmed through the precise localisation of the carbon atoms involved in the 
intermolecular bonds. The spectral features of this phase are revealed through a combination of 
ab-initio lattice dynamics calculations and inelastic neutron scattering experiments. The neutron 
scattering observables are found to be in very good agreement with the simulations which predict 
a partial charge transfer from the Li atoms to the Cgo cage. The absence of a well defined band 
associated with the Li atoms in the experimental spectrum suggests that this species is not ordered 
even at the lowest temperatures. The calculations predict an unstable Li sublattice at a temper¬ 
ature of ~ 200 K, that we relate to the large ionic diffusivity of this system: low frequency optic 
modes of the Li ions couple to the soft structure of the polymer. 

PACS numbers: 61.05.F, 63.20.-e, 63.20.dk, 61.48.-c, 78.70.Nx, 81.05.Tp, 82.35.Lr 
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I. INTRODUCTION 


Solid fullerene polymerization is a well-established phenomenon 1 . Cgo can polymerize 
by undergoing high temperature and high pressure treatments, after ion irradiation or by 
exposure to light 2 . In all these cases, fullerene units invariably connect together by four- 
membered carbon rings, via [2+2] cyclo-addition reactions-! 

The intercalation of light electron donors (he. alkali atoms Li or Na) in the voids of 
the host fullerene lattice can sometimes induce the formation of polymerized structures. In 
this case, the charge transfer from the metal to the highly electronegative Cgo provides the 
necessary chemical pressure for adjacent buckyballs to be close enough together to estab¬ 
lish iutermolecular bonds. A large variety of ID 4 5 and 2D 6 fulleride polymer structures 
was found when intercalated with small-alkali atoms. In addition to the aforementioned 
cycloaddition mechanism, fullcrenes can connect through single C-C bonds. In the lat¬ 
ter category, lithium intercalated fullerides Li x Cgo (x < 6), are remarkable because their 
structure consists of 2D polymeric fullcrenes interconnected by a sequence of single C-C 
bonds and four-membered carbon rings, propagating along two orthogonal directions in the 
polymer plane 40 . Li 4 C 60 belongs to this class of compounds, together with the isostructural 
alkali-earth intercalated Mg 2 C6o polymer, suggesting that the on-ball charge transfer is the 
key factor driving the polymerization arrangement in these systems 7 . The coexistence of 
the two bonding schemes was also observed in the fullerenium salt C6o(AsF 6 ) 2 . In this solid, 
the fullerene units are oxidized to the state Cgg , a rather unstable condition, which results 
in an unusual ID zigzag polymerizatiorP. 

In addition to its novel structure, LijCgo shows a very large ionic conductivity at low tem¬ 
perature, which is an exceptional property for a solid material. The ionic conductivity 
reaches the value of a ~ 1CD 2 S/cm at room temperature, which is comparable to that 
observed in liquid electrolytes, indicating that this material could possibly find applications, 
for example, in fuel cell devices 9 . This physical property originates from the presence of 
intrinsic unoccupied sites interconnected by 3D pathways in the crystalline structure, al¬ 
lowing the diffusion of the Li + ions. The large amplitude movements of the Cgo cages e.g. 
their rotations or radial deformations, can facilitate the Li + ions jumping from one site to 
another. Such dynamical disorder is intrinsically involved in the atomic mechanisms leading 
to oxygen diffusion in certain oxides, like the onset of MoCb free rotations in Bi 2 gMoioOgg 10 . 
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In general the superionic character of a material follows an order-disorder transition which 
unlocks large amplitude movements. In Li 4 C6o, no such transition is observed and the supe¬ 
rionic nature of the material appears progressively in an apparently very ordered fullcrene 
host structure. 

It is possible to break the intermolecular bonds and recover the monomer phase through a 
moderate thermal annealing. The resulting phase shows an unexpected metallic behaviour 11 
which appears to be in contrast with accepted theories indicating the A 4 Ceo compounds (A 
= alkali metal) to be Mott-Jahn Teller insulators 12 13 . 

In the present work we report on the use of neutron scattering to study Li 4 Ceo- The 
structure is investigated by neutron powder diffraction performed at low temperature, i.e. 
in a regime where the Li diffusion is hindered. Rietveld refinement is used to precisely 
localise the carbon atoms of the distorted fullcrene cage -especially those involved in the 
intermolecular bonds- and to reveal the Li positions in the host lattice. The structural 
information is then used to calculate the lattice dynamics of the polymeric system in the 
harmonic approximation. The phonon dispersion curves, the vibrational density of states 
and the thermodynamic properties of this molecular system are discussed. The validity 
of the model is illustrated by the excellent agreement between the neutron scattering ob¬ 
servables derived from simulations and those observed in our inelastic neutron scattering 
experiments (INS). 

Finally, we present the INS results obtained on this system in the monomeric phase at 
high temperature. The experimental data are interpreted with the help of DFT Molecular 
Dynamics simulations (MD-DFT). 


II. EXPERIMENTAL METHODS 

Li 4 C6o samples were prepared by thermal decomposition of 99% isotopically pure 'Li 
azide. Details of the procedure are described elsewhere 10 . Isotopic enrichment was required 
to avoid neutron absorption arising from the natural fraction of 6 Li. Sample handling was 
performed in a controlled atmosphere (high vacuum or Ar glove box with oxygen and mois- 
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ture values below 1 ppm). 

The neutron powder diffraction measurements were performed at the Institut Laue 
Langevin (ILL), using the Super D2B two-axis diffractometer, operating in Debye-Scherrer 
geometry and in high resolution mode (neutron wavelength A = 1.59432 A). A 1.5 gr sample 
was sealed in a cylindrical vanadium can with indium o-ring. Measurements were per¬ 
formed using a Closed Cycle Refrigerator allowing data to be collected at 3.5 K, with an 
accumulation time of 10 hours. 

Inelastic neutron scattering experiments were performed at the ILL, on the thermal time 
of flight spectrometer IN4C, and on the Filter Analyser Spectrometer INIBeF. Additional 
measurements were performed at the Laboratoire Leon Brillouin (CEA Saclay) on the cold 
TOF spectrometer Mibemol, to follow the evolution of the dynamics during the depoly- 
merisation process at high temperature. On IN4C, several incident wavelengths were used 
(1.1, 1.4, 1.5 and 2.4 A), allowing the dynamics to be measured either in Stokes (phonon 
creation) or in anti-Stokes (phonon annhilation) regimes: using short wavelengths, the dy¬ 
namics can be measured at low temperatures in an extended energy range, while using a 
relatively long wavelength, the dynamics can be measured in an extended energy range in 
the anti-Stokes regime, albeit at relatively high temperature. On INIBeF, the measurements 
were conducted at low T (10 K). In both experiments, the 1.5 g sample was sealed inside 
an A1 flat container with an Indium seal. On Mibemol the 1.5 gr sample was sealed inside a 
cylindrical niobium can and heated up to 800 K during 12 hours using a high temperature 
furnace. The incident neutron wavelength was 5 A giving the maximum neutron flux. The 
spectra were recorded as a function of time in order to follow the de-polymerisation process. 


III. SIMULATIONS 

Total energy calculations, geometry relaxations and molecular dynamics simulations were 
performed using the projector-augmented wave (PAW) formalism 14 15 of the Kohn-Sham 
density functional theory®^, w ithin the generalized gradient approximation (GGA), imple¬ 
mented in the Vienna ab initio simulation package (VASP) 18 19 . The GGA was formulated by 
the Perdew-Burke-Ernzerhof (PBE) density functional 20 21 . For total energy calculations, 
the electronic calculations were performed at the gamma point (k=0) and an energy cut-off 
of 499 eV was used. For the geometry relaxations, the energy cut-off was increased by 30% 
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and a 2x2x1 Monkhorst-PackP 2 k grid scheme was used for increased precision. The break 
conditions for the self-consistent field (SCF) and for the ionic relaxation loops were set to 
10~ 5 eV and 10~ 6 eV A -1 , respectively. The latter break condition means that the obtained 
Hellmann-Feynman forces are less than 10~ 6 eV A -1 for the final geometries (futher details 
about the way such calculations are performed can be found in Ref. f23il . Both full (lattice 
constants and atomic positions) and partial (only atomic positions) geometry relaxations 
were performed. Total energies were calculated for the 108 structures resulting from the 
individual displacements (atomic displacement = 0.03 A) of the 18 symmetry inequivalent 
atoms in the monoclinic (I2/m) phase, along the three inequivalent Cartesian directions 
(±x, ±y and ±z). Phonons are extracted from these calculations using the direct method 
as implemented in the Phonon softwar e! 24 * 25 ! The phonon spectra obtained using the fully 
and partially relaxed structures differ only slightly, so that only the results reported in this 
paper concern the fully relaxed structure. 

The charge distribution was calculated from the VASP relaxed geometries using the Bader 
scheme analysis 26 . We used the code bader developed by Henkelmann’s group 27 which gives 
the partial charge on each atom. 

The results of the lattice dynamics calculations will be discussed based on the usual proper¬ 
ties: phonon dispersion curves, total and partial phonon density of states. In particular, a 
weighted phonon dispersion curve image is also discussed: for each phonon \ jq) characterised 
by energy hcjj(q) and momentum hq, an intensity is calculated as: 

= 03(91 ^)ll 2 A) 

with ej(q \ p) being the polarisation vector of the phonon \jq) for atom // in the unit cell. 
This identifies the atoms involved in each mode, and we will call this quantity 7{ /t }(j, q) the 
participation factor of the ensemble of atoms {/i} in the phonon mode | jq). 

Several thermodynamic functions (heat capacity, entropy, free energy and mean square 
atomic displacement) obtained in the framework of the harmonic approximation 28 ' are also 
presented and discussed. These quantities are extracted from the calculations using the 
Phonon Software 24 . 

The neutron observable S(Q, c a) and G(co) are calculated for a powder according to the PALD 
method 23 29 . Based on a set of modes | jq), the coherent dynamical structure factor in the 
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”one-phonon” approximation 30 S co h(Q, ca) is computed for a large number of Q = q + Ghkh 
with random orientation on a dense equidistant ||Q|| grid. Ghki represents the vectors of the 
reciprocal lattice of the crystal. In phonon creation (Stokes): 

S coh (Q,co) oc n(a;,T) + 1 ~ uj(Q)) ( 2 ) 

w j 

with the one phonon form factor: 

_ N d ucoh -W d (Q) _ t 

*jW) = llE - 7W Q- c jWI c 0)II 2 (3) 

<1.1 v"i 

and the Bose thermal population factor n(u>,T ) = (exp(hoj/kgT) — l)^ 1 with fcg the Boltz¬ 
mann constant. Index d refers to the cith atom among the N,d atoms in the unit cell, with 
mass Md and coherent scattering length Wd(Q) is the Debye-Waller factor of atom 

d. In a final step, S co h(Q,uj ) is regrouped according to an adjustable A||<5|| and Acn and 
orientationally averaged to obtain the isotropic function S co h(Q,uj). 

If the incoherent scattering length of carbon is negligible, that of lithium is not and one 
has to account for this scattering in the calculations. The Li total scattering cross sec¬ 
tion is therefore expressed as S(Q,co) = S co h(Q,cu) + S inc (Q,uj ) with the powder-averaged, 
one-phonon, incoherent cross section approximated to: 


A dQ 


,n(cu,T) + l^{b^ c 


y e -2W d {Q) 


cu 


d =1 


Md 


9d{u) 


( 4 ) 


with gdipj) being the partial phonon density of state of atom d: 


9d H = \\ej(q\ d))\\ 2 5(uj - ujj(q)) (5) 

d jg 

and N is the number of unit cells in the crystal (he. 3 NNj is the total number of degrees 
of freedom in the sample). The calculated S(Q,u) spectra were subsequently adapted to 
the experimental conditions (temperature, instrument dependent scattering (Q,oj) range, 
resolution function) to allow for a direct comparison with the data. In particular, the 
INS-weighted spectra were convoluted with a Gaussian with energy-dependent width. The 
calculated and experimental S(Q,u) spectra were therefore treated the same way. In par¬ 
ticular, the grouping and treatment performed in order to derive the generalised density of 
states (GDOS) is identical for both sets of data. 
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The structure of the high temperature monomeric phase of LRCgo was taken from RefJTTl 
As this structure is orientationally disordered, lattice dynamics simulations cannot be per¬ 
formed. We used DFT molecular dynamics simulations to investigate the atomic displace¬ 
ments in this phase. The MD simulations were performed at 800 K using a Nose-Hoover 
thermostat!^ as implemented in the VASP package. The MD time step was set to 1 fs and 
the total simulation time was 3 ps. The neutron weighted density of states were derived 
from the trajectories as the Fourier transform of the atomic velocity auto-correlation func¬ 
tion. This analysis was performed using the nMoldyn package 33 . The simulation box used 
is restricted to the I A/mmm unit cell, i.e. containing only 2 molecular units. As periodic 
boundary conditions were applied, artificial correlations between the rotating Cgo units can 
arise. However these mostly affect the quasielastic scattering region of the neutron spectra, 
which we do not discuss in the present paper. 


IV. RESULTS AND DISCUSSION 

A. Neutron diffraction 

Diffraction data were collected within the angular range 26 = 10° - 145° with a step of 
0.05° (high resolution mode). Rietveld analysis was performed using the GSAS+EXPGUI 
software. The starting peak profile was determined by a preliminary Le Bail pattern decom¬ 
position and the background contribution was modelled with a polynomial curve containing 
21 coefficients. The structural analysis was started without considering the Li ions, due to 
the low (negative) scattering length of 'Li with respect to C. The initial fullcrene arrange¬ 
ment was obtained by considering the carbon positions of the polymeric chains in KCgo 
(where Cgo are connected by [2+2]-cycloaddition bonds), as described elsewhere 31 . Then, 
the fullcrene chains were correctly rotated around their axis and the necessary transfor¬ 
mations were applied to the atomic coordinates, in order to match the correct symmetry 
and the cell dimensions of Li 4 C 60 (a = 9.329(1) A, b = 9.054(1) A, c = 14.984(1)A, (3 = 
90.91(1)°, space group I2/m), with the [2+2] polymerisation running along the 6-axi^ 35 ! At 
this stage, we found that the minimum in the R wp agreement factor was reached by rotat¬ 
ing the fullerene chains by 102° with respect to their arrangement in KCgo- As expected, 
at this stage, the diffraction profile was not well described during the Rietveld refinement 
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(R wp = 7.68%). However, a significant improvement was observed when the positions of the 
carbons C(15), which are involved in the single C-C bonds, and their neighbours (C(5) and 
C(16)) were refined. Although the refinement process was always rather stable, due to the 
high number of observables available (1401), soft constraints were applied to carbon posi¬ 
tions of the whole asymmetric unit and the convergence was reached by slowly decreasing 
the F-factor (penalty factor) at each cycle. 

Further enhancement of the fit (to R wp = 5%) was gained with the refinement of the 
carbon C(l) forming the four-membered carbon ring and its neighbours (C(2) and C(3)). 
Finally, the fractional coordinates of all the remaining ten C atoms of the asymmetric 
unit were refined stepwise. The best agreement between the observed and calculated data 
was obtained for an interfullerene C(15)-C(15) bond of 1.61(2) A. The latter distance is 
significantly smaller than that previously derived from synchrotron data (dc'(i 5 )_c(i 5 ) = 
1.75(2) A®), falling into a more reasonable range if one compares it with a typical C-C bond 
length. Moreover, no anomalous distances were detected in the fullcrene molecule, the bond 
lengths being close to the standard values of 1.54 A and 1.42 A, for sp 3 and sp 2 C-C bonds 
respectively. 

Solid State NMR and DC and AC conductivity 3 ® clearly indicate that the Li ions are 
pinned at special positions in the lattice at the temperature at which the diffraction data 
were collected. This means that no ionic diffusion is expected at 3.5 K. The structural anal¬ 
ysis performed at room temperature using the synchrotron data distinguished two types of 
symmetry inequivalent Li ions, located respectively inside the pseudo-tetrahedral (LR) and 
pseudo-octahedral (Li 0 ) voids of the parent cubic lattice 30 . In the present, low temperature 
neutron investigation, the study of the small alkali ion position was attempted using Fourier 
map analysis, i.e. exploiting the large contrast arising from the negative neutron scatter¬ 
ing length of 'Li. We found that a significant fraction of intercalated alkali ions occupy 
the position (0.44(1) 0 0.76(1)), which corresponds to the small pseudo-tetrahedral voids, 
in agreement with the previous analysis. The insertion of LR atoms with full occupancy 
improved significantly the quality of the fit (R wp = 3.99%). 

In contrast, the Fourier maps did not provide clear evidence of the occupation of the larger 
pseudo-octahedral sites, showing at these positions an unexpectedly poor contrast. The rea¬ 
son for this observation was first ascribed to the combined effect of the small absolute value 
of 'Li scattering length, as compared to C, in addition to the low amount of the alkali ions 



FIG. 1. Powder diffraction pattern of 'L 4 C 60 at 3.5 K. The data are indicated by black crosses, 
the Rietveld fit by the solid red line and the difference by the solid blue line. Ticks marks indicate 
the reflection positions (i? wp = 3.83%, f?F 2 = 1.16%). The fulleride structure is displayed in the 
inset, with the carbon nomenclature. 

in the compound. In fact, if both effects are taken into account, the contribution of Li to 
the scattered intensity turns out to be two orders of magnitude smaller than that of carbon. 
Another hypothesis is that the Lio ions present a certain degree of disorder, e.g. populating 
different octahedral sites in the lattice with similar energy and experiencing slightly different 
local potentials. This hypothesis is supported by DFT calculations-. The spread of Lio ions 
over different positions in the octahedral voids in the lattice is also expected to strongly 
reduce the contrast of the Fourier maps. 

However, a careful analysis of the scattering reveals the existence of negative intensity close 
to the position (0.28, 0, 0.32), which is significantly shifted with respect to the center of 
the octahedral void. The insertion of the Lio alkali ion at this position with a full occu¬ 
pancy further improved the agreement factor ( R wp = 3.83%), suggesting that Lio ions tend 
to move away from the center of the octahedral sites with decreasing temperature. A simi- 
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lar behaviour was recently observed for highly-doped Li fullerene compounds, where Li ions 
were found to form clusters around the center of the tetrahedral positions of the parent cubic 
lattice 37 in the ground-state, leaving the octahedral sites unoccupied. 

In the final configuration, the closest Li-Li distance found is of 2.9(1) A, similar to the atomic 
distance found in Li metal (3.04 A), while the Li + -Ceo contacts reveal a short distance of 
the order of ~ 2.4 A, which is equal to the sum of the ionic radius of Li + (0.7 A) and the 
Van der Waals radius of C (1.7 A). The result of the Rietveld refinement performed on the 
neutron data is displayed in Fig. [lj 


B. Lattice Dynamics Calculations 

1. Geometry optimization and electronic properties 



(a) (b) 


FIG. 2. a) Structure of LRCeo as obtained from geometry optimisation, in the space group I2/m. 
Grey sticks are carbon bonds, magenta and brown spheres are carbon atoms involved in the co¬ 
valent intermolecular bonds: [2+2] bridges and single bonds respectively. Red and yellow spheres 
represent Liy and Lio respectively, b) View along the c axis of one polymeric plane (the Li ions 
have been omitted for clarity).) 
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The geometry relaxation was performed in two steps: the first one consisted of optimizing 
the atomic positions while keeping the lattice parameters at the values obtained by neutron 
diffraction. The experimental structure described above was used as the starting structure. 
The optimized geometry obtained this way was further relaxed by allowing lattice parameters 
to vary as well. The optimized structure was then considered for phonon calculations. We 
have checked that both relaxed geometries gave comparable results in terms of phonon 
dispersion curves and phonon density of states. The full relaxation however gave a better 
result in terms of there being no negative phonon frequencies in any part of the Brillouin 
zone and a slightly better agreement with the neutron derived phonon DOS. The resulting 
structure has lattice parameters a = 9.38 A, b = 9.12 A, c = 15.59 A, with angles a 
= 7 = 90° and monoclinic angle f3 = 91.18°. This gives a AVjV = 5% expansion of 
the lattice with regards to the experimental structure. Relaxed interatomic distances are, 
a cc ( 66) ~ 1.44 A (sp 2 bonds between two hexagons); a cc (65) ~ 1.42 A (sp 2 bonds between one 
hexagon and one pentagon); a cc (sb ) ~ 1.56 A (sp 3 single bond between adjacent molecules), 
a cc([2 + 2]) 1.64 A (sp 3 bond in the [2+2] cycloaddition). In agreement with the diffraction 

investigations, the Li atoms are substantially displaced with regards to the centers of the 
parent tetrahedral and octahedral voids in cubic Cqo. They form atomic planes of atoms in 
the space between two successive polymer sheets (see Fig. [2j) . 


2. Raman and IR modes 

The crystal symmetry point group of LijCgo is C^h- The vibrational representation of 
the system is composed of the non-degenerate A g (R), B g (R), A U (IR) and B U (IR) irreducible 
representations. The phonon modes at point T are either Raman (R) or Infrared (IR) 
active. The complete list of the modes, classified according to their symmetry and frequency, 
are presented in Tables [T] and |TT} The calculated frequencies are compared to those of 
the principal features observed in the Raman, IR and High Resolution Electron Energy 
Loss Spectroscopy (HREELS) spectra reported in the literature 38 41 . Note that HREELS is 
subject to the same selection rules as IR. The association between the calculated and the 
observed modes is based on the best match between their frequency and other arguments 
like, for example, their relative position with regards to common spectral gaps. Compared to 
bulk Cgo, for which only 13 modes are Raman active and 6 are IR active 42 (including lattice 
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modes), resulting in 19 observable modes, the reduced symmetry of the LRCeo lattice, added 
to the deformation of the fullerene cages due to polymerization, increases the number of 
optically active modes to 95 Raman and 94 IR active modes. This renders the interpretation 
of the Raman and IR spectra based only on the frequency match between the experimental 
and calculated features very difficult. A definitive attribution of the observed features would 
require the calculation of the intensity of the Raman/IR active modes which is beyond the 
scope of this paper. 

In this section we restrict the discussion to two types of Raman active modes which are 
of particular interest and which are extensively discussed in the literature. The first one 
concerns the double peak feature located around 120 meV (119 and 121.5 meV in LRCeo) 
which is often associated to the stretching of the [2+2] bond in polymeric Ceo- However 
its observation in the Raman spectrum of the single bonded polymer Na^eo has raised 
a question about this attribution 3940 . In the experimental Raman spectra, this feature 
appears well separated from the preceding Raman peak at lower frequency by a ~15 meV 
gap. In Tab. [TJ these two extreme features are observed at 119 and 102.7 meV respectively. 
Our calculations also predict a 15 meV frequency range populated essentially by IR active 
modes. They create a gap in the Raman spectra between 108 meV and 123 meV. Therefore, 
we attribute the observed 119 meV features to a A g active mode which is represented in 
Fig. [3(c 


This mode shows a complex displacement pattern, composed mostly of single¬ 


bond deformations, but also of [2+2] contributions. In contrast, the modes located at lower 


frequencies imply essentially movements of the atoms forming the [2+2] bridge. In Fig. 3(a) 


we show the B 9 mode calculated at 102.5 meV. Its displacement pattern represents the 
buckling of the [2+2] square-like bridge out of the polymer plane, i.e along the c direction. 
We attribute to this mode the peak observed at 102.7 meV in the Raman spectra. Other 
Raman active modes, implying the deformation of this intermolecular bridge, are calculated 
in this frequency range: a B s mode at 103.8 meV and a A g mode at 104.7 meV, the latter 
mode being associated with the stretching of the bridge along the b direction (see Fig. 3(b)). 


They could also be attributed to the feature observed in the Raman spectra at 102.7 meV. 


In conclusion, our calculations propose that the features in the Raman spectra observed in 
the range around 105 meV are associated to the [2+2] bond deformations, e.g. at frequencies 
significantly smaller than those of polymeric forms of Cgo obtained at high pressure and high 
temperature. Their contribution to the phonon density of states (DOS) is sharp and observed 
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at ~103-105 meV in the total DOS spectrum shown in Fig. [5] The deformation of the [2+2] 
cyclic bond, therefore appears as a highly localised mode featuring a sharp peak in the pDOS 
of the C atoms involved in this bond (see brown line in the top part of Fig. [5]). 

In contrast, the phonon modes involving single bond contributions are calculated to be 
hybridized with other molecule deformations to a large extent. They give a significant 
contribution in the 120-150 meV range. Their contribution to the DOS is relatively smooth 
(see magenta line in the top part of Fig. [5] top). 

The second Raman feature that we discuss in this section concerns the so-called ’’pentag¬ 
onal pinch mode”. It refers to a C 60 mode, of A g symmetry in bulk C 60 , located at 182 meV 
(1469 cnr 1 ) involving tangential displacements of the atoms and in particular, contraction 
of the pentagonal rings of the fullcrene. Its frequency is known to red-shift depending on the 
number of polymeric bonds per molecul e 43 1 44 1 and/or the charge carried by the fullerene 15 
according to the empirical rule: -0.8 meV per excess electron, -0.3 meV per single covalent 
bond and -0.7 meV per [2+2] covalent bond. A mode in the Raman spectra of Li 4 Cgo is 
observed in this frequency range (see Tab.[I| at a frequency of 178.8 meV, i.e. red-shifted by 
~3 meV with regards to the bulk value. Our calculation predicts a A g Raman active mode 


at a frequency of 178.6 meV. This pure cage deformation mode is represented in Fig. 3(d) 


and presents a pattern close to the one expected from the A s (2) ’’pentagonal pinch mode” 
of bulk Cgo- Following the empirical rule described above, we anticipate a charge transfer 
of about 2e per cage 4 ®21, a value substantially lower than that proposed by Macovez et al. 
who derived a value closer to 3e (2.7e) from the analysis of the photoemission spectrum of a 
Li 4 C6o hint 41 . Both predictions however agree on the partial nature of the charge transferred 
from the lithium to the fullcrene cage, a tendency already known in intercalated fullerene 
system^. The analysis of our simulations, using the Bader 2 ®3 approach, indicates a total 
charge of 726.1 electrons per Cgo which means a charge of ~-3e per cage and a partial charge 
transfer from the Li atoms to the Cgo cage which amounts to ~0.75e per alkali ion. This 
value is in very good agreement with the value proposed by Macovez et al. In addition, 
the calculations reproduce well the experimental frequency of the pentagonal pinch mode 
Raman feature. This observation reinforces the difficulty of deriving the excess of charge 
transferred to the fullerene in polymeric form based on the frequency of the Raman feature 
alone. Any use of the above described rule will lead to an underestimate of the on-ball 
charge transfer. 
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141.9 


144.6 

145.7 

146.8 


152.3 

151.6 

157.3 

157.7 

147.6 

148.5 

155.6 


159.9 


160.2 


162.7 

163.5 

165.4 


161.1 


162.0 

161.8 

165.9 


167.6 

167.4 

169.2 


173.0 

172.1 

167.6 


170.1 

169.8 

173.7 


177.9 


184.6 

184.2 

185.7 


175.8 


178.6 

178.8 

184.7 


189.0 

195.6 





184.7 


186.5 

187.2 










TABLE I. Raman active modes frequencies (in meV) calculated for A g (Left) and B y (Right) 
symmetry. The frequencies of the 45 discernible features observed in the experimental Raman 
spectra obtained at room temperature using a red (1.96 eV) laser excitation on a LLCgo polymer 
sample 40 (see also Ref. f39l) are also reported. The experimental frequencies (in italic) are placed 
inside the Exp. columns at positions in the table the closest to those of modes predicted by the 
simulations. 


3. Phonon dispersion curves and Phonon DOS 


Figure [4] shows the dispersion curves of the phonon branches calculated along several 
high symmetry directions in the Brillouin zone. The high symmetry point coordinates are 
in the system of reciprocal vectors (A, B, C) of the crystallographic unit 


given in table 


Iff 
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Calc. Exp. 

A u 

Calc. Exp. 

Calc. Exp. 

Calc. Exp. 

Calc. Exp. 

Calc. Exp. 

18.9 

25.8 26* 

38.0 38* 

11.5 

15.0 

35.1 

39.7 

44.9 

46.2 

38.5 

41.9 

45.2 

47.2 

48.2 0 

57.3 56* 

46.3 

47.5 

48.6 

61.0 

63.7 

65.6 65* t 

49.0 

59.4 

59.6 

68.1 

78.5 

79.8 

66.1 

66.8 

69.6 

80.7 81.5* 

82.6 

83.3 

72.0 71** 

78.4 76* 

79.1 

86.1 

88.4 

89.6 

80.6 

84.1 

86.1 

92.0 

93.4 94* 

94.6 96* 

88.4 

89.0 

90.8 91* 

108.7 

115.8 

118.3 

98.2 99 T 

104.2 

111.0 

120.2 

130.1 

133.2 

115.7 

116.2 

119.3 

143.2 

147.1 14T t 

151.1 151* 

135.1 

137.9 

140.9 

154.2 

157.2 

159.7 

143.2 

148.8 

152.7 

165.2 

166.4 

168.2 

155.3 155.5* 

162.0 

164.2 

171.6 

174.8 

177.2 177* 

164.9 

166.8 

170.1 

185.8 

186.6 

188.9 

174.0 

175.6 

179.1 




181.9 

185.3 

186.6 


TABLE II. Left: IR and HREELS active modes frequencies (in meV) calculated for A u (Left) and 
B u (Right) symmetry. Frequencies of the principal features observed in the experimental IR (from 
Ref. (38j symbol *) and from HREELS spectra (from Ref. 3U symbol ' ) are also reported. The 
experimental frequencies (in italic) are placed inside the Exp. columns at position in the table the 
closest to those of modes predicted by the simulations. 

cell (I2/m). Due to the weak monoclinic character of the unit cell, propagation vectors 
of the form [ 0 , 0 , £] will have major out-of-plane components and refer to phonon modes 
propagating out of the polymer plane. In contrast, propagation vectors of the form [£, 77 , 0] 
refer to phonon modes propagating in the polymer plane. 

The dispersion curve pattern is typical of an anisotropic molecular crystal, its molecular 
nature being revealed by flat ’’molecular” modes in the high frequency range (E >~ 30 meV) 
and the presence of a small 2 meV gap between 26.5 and 28.5 meV. The strong anisotropy 
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(c)A g Raman active mode at 123 meV (992 (d)A ff Raman active mode 
cm -1 ) at 178.6 meV (1441 cm -1 ) 


FIG. 3. Selected Raman Active modes discussed in the text. The Li atoms are omitted for clarity. 


between the in-plane covalent bonds and the out-of-plane van der Waals interactions is 
observed by the rather flat dependence of the molecular inodes along the e.g. T — Z line, 
while a substantially larger dispersion is calculated for e.g. the modes in the [110, 120 
meV] range along the T — Xi direction. These high frequency molecular modes essentially 
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FIG. 4. Left: Dispersion curves along high symmetry directions. Right: Partial dispersion curves 
along high symmetry directions. From left to right, contributions from Li t, Li o and carbon atoms. 
The color scheme represents the value of the participation factor of the atom in the modes: the 
darker the color, the higher the participation (normalized to unity, i.e. a value of one represents 
maximum participation factor). 


Label 

[Qi 5Q21Q3] 

Label 

[Ql,Q 2 ,Q 3 ] 

Label 

[Qi,Q 2 ,Qs] 

Label 

[Qi,Q 2 ,Q 3 ] 

r 

z 

[0,0,0] 

[0,0,1] 

Xi 

V 

[+0,0] 

x 2 

M 

[0,|,0] 
[3.0,i] 

L 

A 

[iio] 

[|A-i] 


TABLE III. High symmetry point coordinates expressed in the system of reciprocal vectors ( A , 
B, C ) of the crystallographic unit cell (12/rn). Note that due to the weak monoclinic (/3 = 91.18°) 
character of the unit cell, the vectors A, B and C are close to being collinear to the single-bonds, 
[2+2] bonds and normal to the polymer plane directions respectively. 

involve carbon displacements, Li ions being largely immobile, as seen from the weak Li 
participation factors in this energy range (see Figj4jLio and Li^) and the Li partial density 
of states (pDOS) (see Fig|5]). As already discussed for other polymeric forms of Cgo, these 
phonon modes share strong similarities with those of bulk Cgo, but the DOS spectrum also 
differs significantly in some places, which can be ascribed to the reduced symmetry of the 
polymeric cage and to the different nature of the hybridization of some C-C bonds at the 
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cage surface®®]. 

The 2D character of the polymeric sheet is reflected into the acoustic region of the dis¬ 
persion curves. In particular, two of the three acoustic branches have a steep dispersion in 
the T — X '2 direction, corresponding to displacements of the cages in the polymer plane. The 
third acoustic branch, on the other hand, shows a much flatter dependence with q close to 
T and presents a significant up-turn for the T — X\ and T — X 2 direction with increasing 
q. This general feature, typical of a 2D system, is well known from graphite 51 and has also 
been calculated for polymeric Xa^Cm) 52 . 

The anisotropy and reduced symmetry of the system is also clearly seen in the lattice 
mode region (E < 15 meV) of the DOS represented in the bottom part of Fig|5j The flat 
contribution of the out-of-plane transverse acoustic mode along the (X 2 VMA) directions 
gives the first van Hove singularity at 5.5 meV. The anisotropy of the polymeric system is 
further reflected by the pseudo gap (4 meV) separating this contribution from the second van 
Hove singularity in the DOS observed at ~ 9 meV, which results from the flat dependence 
of the in-plane acoustic modes in the same region of the Brillouin zone. 

At higher frequencies, a series of flat optical branches is observed in the ~ 10-20 meV range, 
followed by a pseudo gap in the range ~ 22 meV and an intense flat band at ~ 25 meV. Then 
follows the ’’real” gap of 2 meV between 26.5 meV and 28.5 meV featuring no intensity in the 
DOS spectrum. It is well known that the polymeric forms of Cqo have characteristic modes 
that appear in the ’’Gap” region of bulk Ceo ([8-33 meV], see e.g. 49 50 53 and references 
therein). In particular, cage libration modes, with the Cgo cages rotating around their 
principal axis of symmetry, i.e. implying inter-molecular bonds to be slightly deformed, 
are usually found in this range. These modes, located at 3.5 meV in bulk monomeric Cgo, 
are shifted to the 10-30 meV range for polymeric forms of Ceo, depending on the precise 
bonding scheme. These librations are calculated at 25.7 meV (B 9 (4)), 32.3 meV (B 9 (5)) 
and 32.6 meV (A 9 (5)) for LRCeo- Their rotation axes main lie along the c, a and b unit 
cell vectors respectively, as represented in Fig. [6] In monomeric Cq 0 , the first intramolecular 
vibration is a Raman active five fold degenerate II 7 mode at 33 meV. Polymerisation lifts its 
degeneracy and we calculate ’’ghosts” of this mode distributed over a relatively large number 
of modes in the 10-30 meV frequency range, some involving a significant contribution of Li 
displacements, i.e. hybridized with Li modes. An example of such a mode is shown in 
Fig,7(a) which represents a Raman active mode at 16.9 meV (A 9 (3)). Tab. [T] indicates that 
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FIG. 5. Total (black line) and partial density of state of the different atomic species: all C (grey 
line), Li t (red line and area), Lio ( yellow line and area), C(l) -involved in the [2+2] bonds 
- (brown line and area), C(15) -involved in the single intermolecular bonds- (magenta line and 
area). The total density of state was shifted by 0.02 meV -1 for clarity in the low frequency part 
(bottom curves) of the spectrum. The total DOS is not represented in the high frequency range 
(top curves), as it is totally dominated by the C partial DOS in this range. The C(l) and C(15) 
partial DOS are only shown in the high frequency range. 
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a feature in the Raman spectrum of the polymer is observed at 16.5 meV. All the optical 
branches mentioned above have flat dispersion, and they give sharp contributions to the 
DOS. In particular, the flat dispersion of the cage deformation mode in the 17 meV region 
(see Fig j4j right panel) gives a rather intense and sharp feature in the total DOS and C 
pDOS spectra. 



(a)Libration around the c axis at 25.7 (b)Libration around the a axis 


(c)Libration around the b axis at 


meV (B fl (4)) 


at 32.3 meV (B 9 (5)) 32.6 meV (A 9 (5)) 


FIG. 6. The three rotation-like modes and their symmetry. 


The Li phonon modes extend up to a maximum energy of ~ 45 meV. Figjd] and Fig|5] 
show that they are grouped into three energy bands for each Li ion type: centred around 
12, 19 and 36 meV for LR ions, and around 15, 25 and 40 meV for LR ions. These groups 
imply either ferro (in-phase displacement of the two symmetry related ions) or anti-ferro 
translation modes (out-of-phase displacement of the two symmetry related ions) along the 
a, b and c axis respectively (see Fig. | 8 j) . In the lowest frequency group, the modes concern 
both LR and Lio ions, with concomitant orthogonal displacements along the a and c axis, 


as shown on Fig. 8 (a 


The sharp contributions of these modes to the DOS, and the fact that the Li ions do not 
have a significant participation factor in the acoustic bands suggest a quite weak interaction 
between the Li and C atoms, i. e. in good agreement with the partial charge transfer between 
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(a )"H g ghost” mode at 16.9 
meV (A fl (3)) 


(b)”Hybrid mode” at 18 meV (c)”Hybrid mode” at 30.2 meV (A fl (4)) 


(B s (l)) 


FIG. 7. Selected ”H g Ghost” and ’’Hybrid” modes and their symmetry. 


the alkali atoms and the fullerene cages. The large mass difference between the Li and 
the fullerene explains the very weak coupling of the Li ion displacements to the polymer 
vibrations in the acoustic range, with a vanishing participation factor of the alkali atoms 
for non-zero propagation vectors (i.e. except at the L point). A careful analysis of the 
dispersion curves and polarization vectors reveals a certain degree of hybridisation between 
the fullerene modes and the Lip and Li t modes. In particular, it is found that the Li ions 
follow adiabatically the fullerene movements for the libration modes (see Fig. [6]). A coupling 
between the H g derived modes and the Li translations, either in-plane or out-of-plane is 
also observed (see e.g. the modes shown in Fig. [7]). These hybrid modes imply deformations 
of the fullerene cage and are of particular interest as they displace the Li ions far away 
from their equilibrium positions, invariably contributing to the large ionic diffusivity of this 
system. 

A very sharp feature is observed in the Li t pDOS at a frequency of 19.2 meV, which 
originates from the dispersionless bands (see Fig. [I]). No fullerene modes are calculated 
in this region over the entire Brillouin zone. In contrast the Lio related modes in the 25 
meV range show a significant dispersion and a broader feature in the corresponding pDOS, 






(a)Pure Li t at 11 meV (A g (l)) (b)Pure ferro Li-r mode 

(translations along b) at 18.9 


(c)Pure anti-ferro Lio mode 
(translations along b) at 24.3 


meV (A„(l)) 


meV (B s (3)) 


FIG. 8. Selected pure Li modes and their symmetry. 

which originates from the presence of fnllerene modes in this region, e.g. the libration around 
the a axis. 


4- Thermodynamic functions 

As seen in Figj9j the additive combination of low mass and frequency for Li vibrations 
is responsible for their very large mean square displacement (MSD) ( u 2 ), even at relatively 
low temperature. In particular, it is calculated to be much larger than the carbon atom 
equivalent in the temperature range from 0-800 K: (u 2 LiT ) (0 K) = ( u (800 K). The 
alkali sub-lattices appear therefore as soft lattices with large amplitude harmonic vibrations, 
embedded inside a stable carbon matrix, which shows a small degree of flexibility due to 
the soft cage deformation modes. The difference in stiffness and stability of the carbon and 
lithium sublattices is further reflected in the molar heat capacity and entropy represented 
in Figj9| although the fraction of Li atoms in the system is weak, the specific heat at ~ 100 
K originates to a significant extent (close to a third) from Li vibrations, while the classical 
atomic limit of 3 kg is almost reached at 300 K. By contrast, it is still half kg for carbon at 
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400 K. The large amplitude vibrations of Li also imply a large value of their corresponding 
vibrational entropy at relatively low temperatures, (~ 6 cal/mol.K at ~ 200 K, i.e. six times 
larger than that of carbon). The most interesting consequence of this is the Li derived free 
energy which becomes negative at a temperature of ~ 250 K, suggesting that the Li sub¬ 
lattices become unstable above these temperatures, even in the harmonic approximation, 
while the fullerene framework retains its stability. If one includes anharmonic terms, the Li 
atom framework is likely to have melted at such temperatures. As no important change is 
expected for the global structure, the transition from an ordered Li sub-lattice to a liquid-like 
structure is certainly of second order character. 


C. Inelastic neutron scattering 


1. S(Q,lo) maps of the polymer phase 


In Fig. 10 the measured S(Q , lo) of Li 4 C 60 collected at 330 K on IN4 and using an incident 


wavelength of 2.41 A is shown on the top left panel. It is compared to that obtained using 


a 1.4 A incident wavelength at 10 K (right panel). The bottom panel of Fig. 10 corresponds 


to the PALD simulated S(Q,co) in the same conditions (E = hco). The main characteristic 
features of the data are very well reproduced by the simulations, over a wide (Q, u>) range, 
both in term of E and Q dependence of the spectra. In particular, one observes a strongly 
Q-modulated signal, typical of coherent scattering, with maxima of intensity located around 
3.4 and 5.7 A -1 . These features are commonly observed in fullerene systems (see e.g. Ref. I5D 
and reference therein) reflecting the coherent scattering of carbon and high symmetry of the 
fullerene cages. They are typical of modes involving rotations of the cages. At low Q and ui, 
one observes an intense signal emerging form the elastic scattering at 0.75, 1 and 2.2 A" 1 
in the experimental spectra. These features correspond to acoustic branches which extend 
up to ~ 6 meV according to our simulations. It is interesting to observe that a Q-cut of 
the intensity at an energy of 6 meV shows a rotation-like Q dependence with maxima at 
3.4 and 5.7 A -1 . This suggests a strong hybridization between modes of translational and 
rotational character at the Brillouin zone edges, e.g. strong rotation-translation coupling of 
the cages. This is also observed for modes in the intermediate energy range (hoj < 19 meV) 
where both maxima are observed, while the higher frequency modes (twj > 30 meV) have 
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T [K] 


FIG. 9. From top to bottom: atomic mean squared displacement (MSD), molar free energy (F), 
molar vibrational entropy (S) and atomic heat capacity (C„), for all atomic species (black) and for 
C (grey) Li t (red) and Li o ( yellow) atoms. C v is given per atom and in unit of kp, the Boltzmann 
constant. 
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FIG. 10. S(Q,oj ) (E = fiu:) maps of L^Cgo measured at 330 K and 2.4 A (top left panel, with 
negative energy transfer), and at 10 K and 1.4 A (top right panel, with positive energy transfer). 
The corresponding simulated PALD maps are shown in the bottom panel for comparison. The 
spectra were normalized to the integrated intensity in the [16, 19 meV] regions. 

Q dependences with a completely different scheme. 


2. GDOS G(u) in the polymer state 


The Generalised Density of States (GDOS) of Li 4 C6o was extracted from S(Q } co) following 
the standard data treatment described elsewhere 55 ^ using the “incoherent approximation” - 


it is shown in Fig. 12 


Two datasets are shown. One was measured at 320 K using a relatively large incident 
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FIG. 11. Energy cuts (E = hu>) along selected lines in the S(Q,w) map according to Fig. 10: a) 
E = -5 meV, dE = 2 meV; b) E = -9.5 nreV, dE = 4 meV; c) E = 9.5 meV, dE = 4 meV; d) E = 
-17 meV, dE = 4 nreV; e) E = 17 meV, dE = 4 meV; f) E = -25 meV, dE = 6 meV; g) E = 25 
nreV, dE = 6 meV; h) E = -32 nreV, dE = 8 meV; i) E = 32 meV, dE = 8 meV 


wavelength of 2.4 A in anti-Stokes mode. The other was measured at 10 K in Stokes mode, 
which is composed of several datasets, measured on IN4C and INIBeF spectrometers, in 
order to cover the whole spectrum. To match the partial datasets, they were normalized 
in common energy ranges. The 320 K and 10 K spectra show the same characteristics: a 
first feature at 7 meV, followed by a range dominated by a double peak structure at 17 
and 19 meV. In the pseudo-gap region, a feature is observed at 24 meV which seems better 
defined at low temperature, while significantly softened and split at 320 K. This mode was 
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0.04 



FIG. 12. Generalised density of states derived from inelastic neutron scattering data collected at 
320 K (incident neutron wavelength 2.4 A - circles top curve) and 10 K (IN4C for energy lower 
than 60 meV: incident neutrons 2.2 A, 1.4 A and 1.1 A. INIBeF for energy higher than 60 meV 
- squares). Data are compared to calculations (solid lines from top to bottom): GDOS of carbon 
and LR (black line) and total GDOS (blue line). Partial GDOS for C (grey line), Li t (red line 
and area) and Lio ( yellow line and area) are also shown. For clarity, the GDOS were shifted by 
~ 0.005 units, except for the partial GDOS. For the same reason, the calculated spectra are not 
shown for frequency larger than 45 meV, the corresponding spectra in that range being totally 
dominated by the C partial GDOS. 

observed in the HREELS spectrum of a thin him of LLCg <P. At higher frequencies, the 
GDOS features the typical fulleride molecular modes. 

The data are compared to the GDOS calculated for the model developed in the previous 
section. Two spectra are shown (blue and black lines respectively): the first represents the 
total GDOS, i.e. comprising the sum of the C (grey line), LR (red line and area) and Lio ( 
yellow line and area) partial GDOS contributions. The second one represents only the sum 
of the C and LR contributions, which would correspond to the case in which the Lio atom 
vibrations were absent or distributed smoothly over a large spectrum range, i.e. reflecting 
structural disorder. 

The calculated spectra are in very good agreement with the data over the whole frequency 
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FIG. 13. Top: Generalised density of states (GDOS) derived from inelastic neutron scattering data 
collected between 300 K (polymer phase) and 700 K (monomer phase): red = 300 K, yellow = 610 
K, brown = 630 K, grey = 640 K and black = 700 K). Bottom: GDOS extracted from the MD 
simulations at 800 K (see text) in the monomer phase (black solid lines: total GDOS, grey solid 
line: carbon GDOS, red area Lithium GDOS). 


range. In particular, the principal low frequency features of the experimental data are repro¬ 
duced. The first peak of the doublet (at 17 meV) can be attributed to Cgo vibrations, while 
the second peak of the doublet matches a significant feature of the LD partial GDOS. The 
feature at 24 meV can be attributed to a hybrid mode, with carbon and Lio contributions. 
However, the corresponding contribution in the total GDOS is significantly more intense 
than what is observed in the data. The double peak feature is also much less defined in this 
spectrum. A much better agreement with the experimental GDOS is obtained for the black 
spectrum: the 24 meV peak is less intense (with regards e.g. to the intensity at 33 meV) and 
comparable to what is observed, and the double peak feature is much better defined. This 
suggests that the Li 0 contribution to the GDOS is not really observed in the data, reflecting 
the disordered nature of the Li 0 sub-lattice, in good agreement with the diffraction data. 
Polarized neutron investigations would be useful to isolate the smooth, broad distribution 
of Li modes as they allow the coherent (mostly C) and incoherent (essentially Li) scattering 
to be separated. 
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3. Transformation to the monomer phase 


After the Li 4 C 60 sample was heated progressively to 700 K, the polymer phase was trans¬ 


formed into a monomeric phase. Figure 13 shows the evolution of the GDOS with tem¬ 
perature, obtained on the Mibemol spectrometer using a 5 A incident wavelength. It is 
compared to the total and partial GDOS derived from MD simulations. The spectral weight 
is observed to be transferred to low frequency in the monomer phase compared to the poly¬ 
mer phase, in agreement with a change in the intermolecular force constant. In particular, 
the modes in the [10, 30] meV disappear. The simulations show a very good agreement 
with the experimental spectra, in particular concerning the frequency and shape of the first 
intramolecular peak at 33 meV. The GDOS in the gap region (he. in the [10, 30] meV 
range) is however not nul. According to the simulation, it contains contributions from the 
Li atoms, albeit more structured in the data than in the simulations. Phase segregation into 
Li and Cgo monomer clusters is probably responsible of this excess structured features in 
the gap. 


V. CONCLUSIONS 

In this paper we have reported an extensive neutron scattering investigation of the LRCgo 
compound. Diffraction data defines a structure with a certain degree of disorder associated 
with one kind of Li atom, referred to as Lio 

The lattice dynamics simulations give phonon dispersion curves and phonon DOS and iden¬ 
tify Raman and IR active modes. The PALD method allows the S(Q,oj) maps of the system 
and the Generalysed Density of States to be generated and compared with neutron scat¬ 
tering observables - the agreement is very good. The DFT simulations are coherent with 
a 3e charge transfer from the Li to the Ceo molecule and the peculiar polymeric bonding 
scheme proposed initially from diffraction and Raman investigations. The inelastic features 
from the Lio atoms are not observed in the experimental data, suggesting a disordered Lio 
sublattice. The presence of hybrid Li-C modes at relatively low frequencies, associating Li 
translations and either Cgo rotations or cage deformations, are highlighted. These modes 
have to be considered as part of the atomic description of the ionic diffusivity in this solid. 
In particular, a peak in the GDOS at 24 meV, attributed to an hybrid C and Lio feature, 
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is found to be temperature dependent. The detailed temperature dependence of this peak 
should be used as a probe of the influence of the Li disorder and large amplitude C 60 modes 
on the ionic conduction in this system as it could reveal the microscopic origin of the Li 
mobility of this crystal. 
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